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Abstract 


The zonal mean atmospheric cloud radiative effect, defined as the difference of the top- 
of- atmosphere (TOA) and surface cloud radiative effects, is estimated from three years of 
Clouds and the Earth’s Radiant Energy System (CERES) data. The zonal mean short- 
wave effect is small, though it tends to be positive (warming). This indicates that clouds 
increase shortwave absorption in the atmosphere, especially in midlatitudes. The zonal 
mean atmospheric cloud radiative effect is, however, dominated by the longwave effect. 
The zonal mean longwave effect is positive in the tropics and decreases with latitude to 
negative values (cooling) in polar regions. The meridional gradient of cloud effect between 
midlatitude and polar regions exists even when uncertainties in the cloud effect on the 
surface enthalpy flux and in the modeled irradiances are taken into account. This indi- 
cates that clouds increase the rate of generation of mean zonal available potential energy. 
Because the atmospheric cooling effect in polar regions is predominately caused by low- 
level clouds, which tend to be stationary, we postulate that the meridional and vertical 
gradients of cloud effect increase the rate of meridional energy transport by dynamics in 
the atmosphere from midlatitude to polar region, especially in fall and winter. Clouds 
then warm the surface in polar regions except in the Arctic in summer. Clouds, therefore, 
contribute in increasing the rate of meridional energy transport from midlatitude to polar 
regions through the atmosphere. 
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1. Introduction 


Clouds perturb the net top-of-atmosphere (TOA) irradiance from that of clear-sky condi- 
tions (e.g. Ramanathan et al. 1989). The spatial distribution of radiation deposited to 
the earth generates the temperature gradient. General circulations reduce the meridional 
temperature gradient by transporting the energy poleward. Because dynamics is driven by 
the temperature gradient, it is postulated that poleward energy transport by dynamics is 
likely to be affected by the presence of clouds. Stuhlmann and Smith (1988) showed, using 
climatological low and mid level clouds, that both cloud types contribute to an increase 
in the rate of generation of zonal available potential energy in mid-latitude and polar re- 
gions. Because only a part of available potential energy is converted to kinetic energy 
(Peixoto and Oort 1992), this does not necessarily means that clouds increase the rate of 
poleward energy transport. Zhang and Rossow (1997) estimated the effect of clouds on 
meridional energy transport from the International Satellite Cloud Climatology Project 
(ISCCP, Schiffer and Rossow 1983) global radiative flux dataset. They concluded that 
clouds enhance meridional energy transport by the atmosphere and reduce the transport 
by the ocean. Zhang and Rossow (1997), however, did not consider the surface latent heat 
and sensible heat flux effects in estimating cloud effects on meridional energy transport. 
Solin and Smith (1992) discussed the effect of cloud type on the meridional circulation 
but they did not. explicitly compute the zonal mean atmospheric cloud effect. Randall et 
al. (1989) used a general circulation model to estimate the effect of the atmospheric cloud 
radiative effect. They concluded that the atmospheric cloud radiative effect over oceans 
intensify Hadley circulation, tropical easterlies, and subtropical westerly jets. 
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In this study, we use TOA and surface shortwave and longwave irradiances from Clouds 
and the Earth’s Radiant Energy System (CERES, Wielicki et al. 1996) data and estimate 
the zonal mean cloud radiative effect at the TOA and surface and to the atmosphere. We 
also include the cloud effect on the surface latent and sensible heat fluxes (enthalpy flux) 
using the National Center for Environmental Prediction (NCEP) reanalysis (Kalnay et al. 
1996) to analyze zonal mean atmospheric cloud effect on the meridional energy transport by 
the atmosphere quantitatively. We then investigate the atmospheric cloud effect by cloud 
height and optical thickness so that the effect on meridional energy transport, especially 
from midlatitude to polar regions, can be understood. 

2. Method 

2.1. CERES Data and Daily and Monthly Averaging 

CERES data from the Terra satellite taken from March 2000 through February 2003 are 
used in this study. TOA and surface irradiances are computed by a two-stream model (Fu 
and Liou 1993; Kato et al. 2005) by forcing the agreement of the TOA irradiance with 
that derived from CERES radiance measurements by angular distribution models (Loeb et 
al. 2005; Kato and Loeb 2005). Inputs to the two-stream model are MODIS derived cloud 
and aerosol properties, MODIS-derived surface skin temperature, ozone amounts (Yang 
et al. 2000), and temperature and relative humidity profile from the Goddard EOS Data 
Assimilation System (GEOS-4, Bloom et al. 2004). Six hour and 1° x 1° GEOS-4 maps 
are linearly interpolated in space and time for the CERES overpass time and location. 
When MODIS-derived aerosol properties are not available, modeled aerosol properties by 
a transport model (Model of Atmospheric Transport and Chemistry, MATCH Collins et 
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al. 2001) are used. The ocean spectral surface albedo is from Jin and Stamnes (1994). 
Broadband land surface albedos are inferred from the clear-sky TOA albedo derived from 
CERES measurements (Charlock et al. 2006). More detailed descriptions of the irradiance 
computation in the CERES project are found in Charlock et al. (2006). TOA modeled 
irradiances are constrained to match CERES-derived irradiances by altering mainly cloud 
properties by the method described in Rose et al. (1997). Those instantaneous irradiances 
are included in a CERES product called Cloud and Radiative Swath (CRS). The edition 
2B, which is available from the NASA Langley Atmospheric Science Data Center, is used 
in this study. 

We convert instantaneous irradiances to daily mean irradiances by the method de- 
scribed in Appendix A. The underlying assumption in the process is that meteorological 
conditions do not change over the course of a day. The error in the daily value, there- 
fore, depends on the actual diurnal variation of albedo and transmittance. For example, 
a systematic diurnal variation of water vapor or cloud amount that is not sampled at the 
overpass time (~ 10:30 AM) is not considered in our estimate. In addition, the solar zenith 
dependent transmittance is estimated by averaging instantaneous transmittance computed 
for the overpass time as a function of solar zenith angle and scene type. Because Terra 
is on a sun-synchronous orbit, the transmittance for large solar zenith angles is from high 
latitude regions. Surface irradiance near sunrise and sunset at low latitude regions is, 
therefore, computed using transmittance derived from high latitude regions. As a conse- 
quence, the daily downward surface irradiance in low latitude regions is overestimated. In 
addition, no effort is made to correct the instantaneous longwave irradiance for the daily 
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value; we simply averaged daytime and nighttime values weighted by the day-night frac- 
tion. We estimate the error in shortwave and longwave surface irradiances by comparing 
with surface observations in sections 3.2 and 3.5. 


3. Results 

3.1. Zonal Net Irradiance 

We define the net shortwave and longwave irradiances at the top of the atmosphere and 
surface as the downward minus upward irradiance, 


rpx _ t?x 
sw sw.dn 


_ J? x 

sw,upi 


( 1 ) 


and 


rpx _ rpx -px 

r Iw r lw,dn J -lw,up'> 


( 2 ) 


where the superscript x is either sfc or TO A. The net irradiances are, therefore, positive 
when the radiative energy is deposited to the system. The net irradiance F ra d is the sum 
of the net shortwave and net longwave irradiances, 


rp x _ rpx I rx /o\ 

r rad r sw ' r Iw 

Net radiative energy deposition to the atmosphere is given by the net irradiance at TOA 
minus the net irradiance at the surface, 

p atm _ pTOA _ psfc (4) 

where the subscript y is sw, Iw, or rad. 

The net shortwave irradiance of the atmosphere , which is the shortwave ab- 
sorption by the atmosphere, is primarily a function of solar zenith angle (Figure la). The 
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seasonal mean net shortwave irradiance of the atmosphere shown in Figure la increases 
toward the tropics. However, the absorbed irradiance by the atmosphere in the summer 
hemisphere is nearly constant with latitude and rarely exceeds 100Wm~ 2 (Figure la). 
When is divided by the daily mean insolation at TOA, the seasonal mean shortwave 

absorptance of the atmosphere is between 0.2 and 0.25 except in polar regions (Figure 
Id). When the absorbed irradiance is averaged, the global and annual mean shortwave 
absorptance of the atmosphere is 0.217. This is significantly larger than the absorptance 
of 0.196 estimated by Kiehl and Trenberth (1997) and the absorptance of 0.208 estimated 
by Raschke et al. (2005) from ISCCP-FD data (from 1991 through 1995). As mentioned 
above, our estimate is based on a two-stream model using cloud properties derived from 
MODIS and water vapor profile from GEOS 4 while Kiehl and Trenberth (1997) used 
a standard atmosphere with 49% low-level, 5%, middle-level, and 20% high-level clouds. 
Aerosols are included in our estimate but excluded in the estimate by Kiehl and Trenberth 
(1997). The difference is then attributed to the difference in cloud properties (optical 
thickness, height, droplet size, phase), aerosol properties, and water vapor amount. The 
net longwave irradiance of the atmosphere Fj^ m is more negative in the tropics than in 
polar regions (Figure lb). Fj^ m shows a larger seasonal variation in the northern hemi- 
sphere than the southern hemisphere (Figure lb). Because the larger seasonal variation in 
the northern hemisphere also appears in the clear-sky net longwave TOA irradiance (Fig- 
ure le), it is likely caused by the seasonal variation in the water vapor amount. The net 
longwave irradiance F^ m over the Antarctic is significantly less negative than Fj^ m over 
the Arctic because of the lower temperature and water vapor amount over the Antarctica. 
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When the net shortwave and longwave atmospheric irradiances are added, F^*™ is 
negative (Figure lc) because the magnitude of Ff^™ is larger than the magnitude of F^ m . 
The value is about -lOOWm -2 for all latitude but F“*™ is more negative in the southern 
hemisphere than in the northern hemisphere. Our estimate of the annual mean atmospheric 
net irradiance is similar to that given by Zhang and Rossow (1997) from ISCCP FC data, 
although their estimate does not show the difference between the northern and southern 
hemispheres. 

3.2. Comparison with Surface Observations 

As mentioned earlier, our estimate of surface irradiance is based on a radiative trans- 
fer model. Assumptions made in the estimate process introduce errors in the daily and 
monthly mean shortwave and longwave surface irradiances. In this section, we compare 
modeled surface irradiances with observations for the error estimate. Note that the mod- 
eled shortwave and longwave TOA irradiances are tuned to agree with CERES irradiances. 
As a result, the daily global annual mean TOA irradiance agrees with that derived from 
CERES measurements to within 0.7% for shortwave and to within 0.2% for longwave. 

We use observations taken at three sites, Manus (2°S 147°E, tropics), Southern Great 
Plains (36°N 97° W, mid-latitude), and Barrow (71°N 156° W, Arctic). Radiation data 
at Manus and Southern Great Plains were taken as a part of the Atmospheric Radiation 
Measurement (ARM, Stokes and Schwartz 1994; Ackerman and Stokes 2003) program run 
by the Department of Energy. The NOAA/CMDL Solar and Thermal Radiation (STAR) 
group took radiation data at Barrow, Alaska. Observed shortwave and longwave irradi- 
ances taken from March 2000 through Feb. 2003 are averaged to compute monthly mean 
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values. Figure 2 shows a comparison of monthly mean shortwave and longwave surface 
downward ir radiances and Table 1 lists the annual mean shortwave and longwave irradi- 
ances for the three sites. The modeled annual mean downward shortwave irradiance over 
three 1° x 1° grid boxes which contain the three observation sites are 7.8% greater than 
the corresponding observations when the differences from the three sites are averaged. 
Similarly, the mean relative difference of the modeled annual mean downward longwave 
irradiance over three 1° x 1° grid boxes containing the three observation sites and cor- 
responding observed irradiances is 1.1%, where modeled irradiances at all three sites are 
smaller. When modeled shortwave and longwave irradiances are compared with obser- 
vations at 25 Baseline Surface Radiation Network (BSRN, Ohmura et al. 1988) sites, 
including above three sites, the mean relative difference is 6.5% for shortwave and -3.2% 
for longwave and the standard deviation is 5.8% and 5.1%, respectively. Because the dif- 
ference from the three sites is within the standard deviation of all BSRN sites, the above 
three sites are representative for all BSRN sites. Therefore, we use values from the three 
sites to estimate the errors in the modeled atmospheric irradiance below and Section 3.5. 
While the reason for the difference of modeled and observed downward surface irradiances 
especially for shortwave needs to be investigated in future, we treat the difference as the 
error in the model in this study. 

Upward shortwave and longwave irradiances are more problematic because measure- 
ments over a small area do not represent a larger area that CERES instruments observe. 
This is especially true for the Manus site since the most of area within the 1° x 1° grid 
box is ocean. The mean relative difference of the modeled and observed upward shortwave 
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irradiance of SGP and Barrow sites is 1.7% (modeled irradiance is greater) and the mean 
relative difference of the upward longwave irradiances of the two sites is 0.9% (modeled 
irradiance is greater). These differences lead to 16.1% under estimate net shortwave irradi- 
ance of the atmosphere and 4.4% overestimate net longwave irradiance of the atmosphere 
when the differences over the three sites are averaged. The relative difference of short- 
wave plus longwave net irradiance of the atmosphere is -10.8%. Table 1 summarizes the 
difference between modeled and observed annual mean irradiances at three sites. 

3.3. Comparison with IceSat-derived Cloud Cover 

Because the surface irradiance is computed with MODIS derived cloud properties, it is 
affected by retrieval errors such as those in cloud fraction, optical thickness, cloud height, 
droplet size, and phase. A part of error in the surface irradiance shown in the previous 
section is, therefore, caused by the error in retrieved cloud properties. Even though com- 
parisons of these properties with those derived from more accurate methods are required 
to fully understand the cause of the error, we expect that the surface irradiance is most 
sensitive to the cloud fraction among cloud properties. We, therefore, only compare the 
cloud fraction with that derived from more accurate measurements in this section. 

Figure 3 shows the meridional cloud fraction derived from MODIS (Minnis et al. 2003) 
and from the Geoscience Laser Altimeter System (GLAS) on IceSat (Zwally et al. 2002). 
The GLAS cloud fraction is derived from a 1064nm laser (medium resolution, GLA09 re- 
lease 26) and cloud optical thickness is derived from a 532 nm laser (medium resolution, 
GLAll release 24). The maximum optical thickness that GLAS can derive is 5. The 
monthly mean cloud fraction over northern hemisphere midlatitude and southern liemi- 



sphere midlatitude derived from GLAS is 0.73 and 0.78, respectively. The corresponding 
cloud fraction from MODIS derived by the CERES cloud algorithm is 0.58 and 0.79, re- 
spectively. The threshold of cloud detection of the CERES cloud algorithm is 0.3 (Minnis 
et al. 2003). The fraction of clouds having optical thickness greater than 0.3 derived from 
GLAS is 0.69 for the northern hemisphere midlatitude and 0.75 for the southern hemi- 
sphere midlatitude. If we neglect the sampling time difference between two satellites, the 
error in the MODIS derived cloud fraction is less than 0.1 in most regions when the cloud 
fraction from GLAS-derived optical thickness greater than 0.3 is compared. 

3.4. Zonal Mean Cloud Effect 

To understand cloud effects on the energy budget of the atmosphere, we compute the zonal 
mean cloud radiative effect at the top of the atmosphere and surface, as well as the zonal 
mean cloud radiative effect to the atmosphere. The cloud radiative effects are computed as 
the net irradiance under all-sky conditions minus the net irradiance under clear-sky condi- 
tions (Ramanathan 1987; Ramanathan et al. 1989). A positive value indicates a warming 
effect by clouds. To avoid the influence of the water vapor amount and aerosol property 
differences between clear-sky and all-sky conditions (Li and Trishchenko, 2001), we use 
computed irradiances with and without clouds for the estimate. In addition, this approach 
provides the cloud radiative effect for all samples as opposed to the approach using the dif- 
ference between observed all-sky and clear-sky irradiances. Note that both shortwave and 
longwave irradiances computed with and without clouds are included in the CRS CERES 
product. The zonal mean shortwave effects at TOA and the surface are both negative (Fig- 
ure 4). The zonal mean shortwave effect to the atmosphere is generally positive, which 
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indicates that the zonal mean cloud radiative effect at the surface is more negative than 
that at TOA. A large radiative effect to the atmosphere occurs in midlatitudes in spring 
and summer in both hemispheres. The zonal mean atmospheric cloud radiative effect is 
slightly negative (less than 1 Wm~ 2 ) in the tropics and polar regions (10Wm“ 2 ). Note 
that the negative zonal mean shortwave cloud effect in tropics is within the error estimated 
in the previous section (the model underestimates by ~ 10Wm -2 ), if all bias error occurs 
under cloudy conditions. The zonal mean atmospheric shortwave cloud effect is, therefore, 
mostly positive. Low-level clouds increase the shortwave absorption in the atmosphere by 
increasing the photon path length while high-level clouds reduce the transmittance to ab- 
sorbing water vapor layers at a lower part of the atmosphere. Figure 4 indicates, therefore, 
increasing absorption by low-level clouds is larger than decreasing absorption by high-level 
clouds when the effects are weighted by the corresponding probability of cloud occurrence. 

The zonal mean cloud longwave effect is shown in Figure 5. When viewed from TOA, 
clouds reduce the effective temperature by reducing the atmospheric transmittance and 
increasing emission from the altitude where clouds are located. The effective temperature 
difference from the clear-sky value decreases as the surface temperature and cloud height 
decrease so that the zonal mean TOA cloud longwave effect decreases with latitude (Figure 
5a). A smaller zonal mean TOA longwave cloud effect near 20N and 20S is due to a small 
amount of clouds in these regions. When viewed from the surface, clouds increase the 
emissivity of the atmosphere and effective temperature. The emissivity difference from a 
clear-sky value increase as water vapor amount decreases so that the zonal mean cloud 
radiative effect at the surface increases with latitude (Figure 5b). The drop in the zonal 
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mean surface cloud effect over Antarctica is also due to a smaller cloud amount compared 
to the cloud amount near 60S. When the zonal mean surface effect is subtracted from 
the zonal mean TOA effect, the zonal mean atmospheric longwave effect is positive in 
the tropics and almost linearly decreases with latitude to negative values in polar regions 
(Figure 5c). This simple linear relation with latitude is caused by the combination of 
the latitudinal dependence of the difference between the mean cloud top height and water 
vapor effective emission height determined from the surface and the latitudinal dependence 
of the cloud fraction (Appendix B). 

When zonal mean shortwave and longwave cloud radiative effects are combined, clouds 
reduce the net irradiance at TOA except in midlatitudes in the winter hemisphere (Fig- 
ure 6a). The zonal mean surface cloud effect is negative in the tropics and midlatitude 
and positive in polar regions except in the Arctic summer (Figure 6b). The shortwave 
effect, therefore, dominates in the cloud radiative effect at TOA and surface in the sum- 
mer hemisphere. Because the shortwave effect at the TOA and at the surface are both 
cooling and nearly equal, the longwave effect dominates in the atmosphere (i.e. Figure 6c 
closely resembles Figure 5c). Subsequently, when the zonal mean shortwave and longwave 
effects are combined, a gradient in the meridional atmospheric cloud radiative effect exists, 
positive in the tropics and negative in polar regions. The meridional gradient persist for 
all seasons. Because the sign changes between low and high latitudes, the global mean 
radiative cloud radiative effect to the atmosphere is only 2.2 Wm~ 2 when the zonal mean 
atmospheric cloud radiative effect is averaged. Although the sign of the cloud radiative 
effect is different, this small annual and global mean net atmospheric cloud radiative effect 
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is consistent with the result by Raschkc et al. (2005) who show that the global mean value 
from ISCCP data is -1.3 Win -2 . 

The cloud radiative effect depends on cloud properties. To understand the atmo- 
spheric cloud radiative effect dependence on cloud properties, Figure 7 shows the net 
cloud radiative effect as a function of cloud optical thickness and cloud top pressure. Only 
single-layer clouds are included in Figure 7 to show the dependence clearly. Cloud optical 
thickness primarily affects the shortwave irradiance while cloud height primarily affects 
the longwave irradiance. Note that the atmospheric cloud radiative effect shown in Figure 
7 (second row) is the effect in the entire atmospheric column as a function of cloud top 
height and cloud optical thickness; it is not the vertical profile of the effect. Figure 7 shows 
that the atmospheric cooling effect in polar regions is caused by low-level clouds, which 
occur most frequently in polar regions (Figure 7 bottom row). The atmospheric effect of 
high-level clouds is warming for all regions (Figure 7). The warming effect is larger in the 
tropics than in polar regions because the frequency of occurrence of thick high-level clouds 
is larger and the mean cloud top height of high-level clouds is higher in the tropics. 

In addition to radiation, clouds also expect to affect the enthalpy flux (i.e. sensible 
heat and latent heat fluxes) at the surface. Because our goal is to estimate the cloud effect 
on meridional energy transport by the atmosphere, we need to quantify the cloud effect 
on the surface enthalpy flux. We use the surface enthalpy flux from the NCEP reanalysis 
(Kalnay et al. 1996) for the estimate. We assume that the enthalpy flux is proportional 
to the net surface irradiance so that the cloud effect on the enthalpy flux Fh — Fn, c ir is 

F„ - F„' C ir = F„( 1 - Flll'JF’ll), (5) 
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where Fjj is the enthalpy flux under all-sky conditions, F is the net irradiance at the 
surface and the subscript clr indicates the clear-sky condition. The cloud effect F*jJ . 
combined the radiative and surface enthalpy fluxes, is, therefore, 


J? s f c _ I ? s f c _ T? s f c — j T? s f c 
r cld ~ r r clr ~ \ r rad 


+ F S H fC )(Fi 


rad 


_ p s f c )/p s f c 

rad, clr ) / rad * 


(6) 


With this assumption, the total cloud effect is the cloud radiative effect multiplied by the 
ratio of the sum of the net irradiance and enthalpy flux to the net irradiance. 

The parameterization of surface enthalpy flux suggests that the flux is a function 
of wind speed and local gradient between sea surface and atmosphere (Washington and 
Parkinson 1986, pl22, Fairall et al. 1996). We need to test, therefore, if there is any relation 
between the net surface irradiance and surface enthalpy flux using observations. For this 
purpose, we used the surface enthalpy flux from Hamburg Ocean Atmosphere parameters 
and fluxes from Satellite monthly data (HOAPS-3, 0.5° global grid, Grassl et al. 2000). 
Figure 8 shows the surface enthalpy flux anomalies divided by the monthly mean value 
as a function of the net surface irradiance anomalies divided by the monthly mean value 
from March 2000 between 20° N and 20° S. The anomalies are defined as the deviation 
of the zonal mean values from the monthly mean value over the four-year period from 
March 2000 through Dec. 2004. Although, Figure 8 indicates that the surface enthalpy 
flux change is positively correlated with the surface net irradiance change, the correlation 
coefficient is 0.16. If ocean surface property changes are caused by radiation and surface 
flux is balanced by the surface net irradiance, the surface enthalpy flux is likely to correlate 
with the surface net irradiance. A weak correlation shown in Figure 8 indicates, however, 
that the surface enthalpy flux anomalies are not predominately driven by the net surface 
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irradiance anomalies over tropical oceans in a monthly time scale. We, therefore, consider 
that the cloud effect on the surface enthalpy flux estimated with (5) is the maximum cloud 
effect on the surface enthalpy flux. The minimum effect is simply no cloud effect on the 
surface enthalpy flux. The uncertainty envelope of the atmospheric cloud effect is then 
given with and without the assumption (5) on the surface enthalpy flux. The uncertainty 
in the cloud effect on the surface enthalpy flux is large in the tropics where the surface 
enthalpy flux is large. When the surface enthalpy flux is small, such as in midlatitudes 
and polar regions, the uncertainty envelope is also small. 

When the cloud effect on the surface enthalpy flux is included with the assumption 
that the cloud effect on the enthalpy flux is proportional to the net surface irradiance, the 
positive effect in the tropical atmosphere is almost eliminated (Figure 9). Therefore, the 
uncertainty in the cloud effect on the surface enthalpy flux is too large to determine the 
atmospheric cloud effect observationally over tropics. However, the meridional gradient 
of atmospheric cloud effects between mid-latitude and polar regions remains because of a 
stronger cooling caused by clouds in polar regions and a smaller net enthalpy flux at the 
surface. 

Our zonal mean atmospheric cloud radiative effect in the tropics is smaller than that 
estimated by Randall et al. (1989, Figure 20) who ran a general circulation model with 
and without clouds over oceans to estimate zonal mean atmospheric cloud radiative effect. 
Their result indicates that the cloud radiative effect is to intensify Hadley circulation 
and tropical easterlies. These are caused by a larger tropical boundary layer wind speed, 
stronger surface evaporation, and a larger water vapor amount in the atmosphere. In their 
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study, therefore, the cloud effect on the surface enthalpy flux is to increase the flux, further 
warm the atmosphere and cool the surface instead of offsetting the cloud radiative effect 
in the tropics. This is opposite in effect compared with the cloud effect estimated with 
our assumption. Therefore, their zonal mean atmospheric cloud effect is larger than our 
estimate in the tropics. This suggests that the cloud effect on the surface enthalpy flux by 
(5) is too large and the gradient of zonal mean cloud effect between the tropics and mid- 
latitude might exist. Our zonal mean atmospheric cloud radiative effect over midlatitudes 
and polar regions is, however, similar to that estimated by Randall et al. Our result- 
of the meridional gradient of cloud effect in mid-latitude and polar regions is, therefore, 
consistent with the result by Randall et al. (1989) that clouds intensify subtropical westerly 
jets because the meridional gradient of the atmospheric cloud effect between midlatitude 
and polar regions can enhance the vertical shear of geostrophic winds through the thermal 
wind relation (e.g. Holton p75). 

3.5. Meridional Gradient Error Estimate 

If we take the mean net irradiance of the atmosphere from three sites used for the compar- 
ison, the relative difference of -10.8% estimated in Section 3.2 corresponds to a 9.5Wm~ 2 
bias error. As shown in Table 1, the net irradiance of the atmosphere over the SGP site is 
biased high by 3.1 Win -2 and that over the Barrow site is biased low by 9.7 Wm -2 . If we 
assume that all the error occurs in cloudy conditions and the error in clear-sky irradiances 
are negligible, the bias error in the zonal mean net atmospheric radiative cloud effect is 
-9.5Wm~ 2 and the uncertainty in the gradient between midlatitude and polar regions is 
12.8 Wm~ 2 . This error estimate in the zonal mean cloud effect is the upper limit because 
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the error in the shortwave daily mean irradiance is partly due to underestimate water vapor 
amount in the atmosphere at large solar zenith angle in low latitudes, which affects both 
all-sky and clear-sky downward shortwave irradiances. A comparison with IceSat derived 
cloud fraction shows that the MODIS-derived cloud fraction is biased low at northern 
hemisphere midlatitude and biased high at 70°N. Because the atmospheric cloud effect in 
polar regions is more cooling than that in midlatitudes, the cloud fraction error at two 
sites is consistent with the error estimate of the meridional gradient of cloud effect. Note 
that the error estimate derived from a comparison with surface observation includes the 
effect of errors in cloud retrieval since modeled surface irradiances were computed with 
MODIS-derived cloud properties. 

Zhang and Rossow (1997) estimated that the uncertainty in the zonal mean surface 
enthalpy flux is at least 20Wm 2 . This corresponds to about 13% of the net surface irradi- 
ance over tropical oceans. This leads to a 13% uncertainty in the surface cloud effect by (6), 
which gives about a 5Wm -2 uncertainty in both the zonal mean surface and atmospheric 
cloud effect. Because of a smaller surface enthalpy flux in midlatitude and polar regions 
than in tropics, we expect that the uncertainty in the zonal mean atmospheric cloud effect 
due to the uncertainty in the surface enthalpy flux is smaller than 5Wm~ 2 in midlatitude 
and polar regions. Therefore, the gradient of the cloud effect between midlatitudes and 
polar regions shown in Figure 9 is significant even when the modeling error and uncertainty 
in the surface enthalpy flux are considered. More importantly, the meridional gradient of 
cloud effect is physically plausible as described in section 3.2. We therefore conclude that 
the meridional gradient of cloud effect in midlatitude and polar regions is significant and 
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subsequently estimate the cloud effect on meridional energy transport in section 4. 


4. Discussion 

Clouds present in polar regions cool the atmosphere more than those in mid latitude and 
tropics cool their residing atmosphere. The meridional gradient of the atmospheric cloud 
effect enhances the temperature gradient in the atmosphere from that found under clear- 
sky conditions. Because mean zonal available potential energy is generated by heating 
warm air at low latitudes and cooling cold air at high latitudes (Lorenz, 1955; Peixoto and 
Oort 1992, p377), this indicates that clouds increase the rate of generation of mean zonal 
available potential energy in the atmosphere (Stuhlmann and Smith 1988). 

Energy transport to polar regions by dynamics can compensate for cooling by clouds 
if zonal mean available potential energy is converted to kinetic energy. If energy transport 
by dynamics in the atmosphere does not compensate cooling by clouds, the atmosphere is 
simply colder than the atmosphere with no clouds and the meridional temperature gradi- 
ent in the atmosphere can intensify with time without changing the TOA net irradiance 
significantly. Because clouds alter the net TOA irradiance as shown in Figure 6, it is pos- 
tulated that clouds affect the rate of meridional energy transport (e.g. Zhang and Rossow 
1997; Moore and Vonder Haar 2001; Weaver 2003). In this paper, we therefore, estimate 
the rate of energy transport by the atmosphere equivalent to the zonal mean cloud effect 
using a CERES data set. We then examine, based on the zonal mean atmospheric effect 
by cloud type, whether or not dynamics can indeed utilize the meridional gradient of the 
cloud effect in transporting the energy poleward. 

Before a detailed analysis, the significance of the meridional gradient of the zonal 


17 



mean cloud effect can be understood from a simple scale analysis. If we take 


A TT'atm 
^ old 

Lc p ph 


( 7 ) 


as the meridional temperature gradient generated by clouds per unit time, where AF^J 71 
is the difference in the cloud effect separated by the horizontal length L, p is the density of 
air, h is the depth of the atmosphere, and c p is the specific heat capacity of air at constant- 
pressure, and use AF ra( j — 10Wm -2 over 1000 km, (7) is approximately lx 10 _12 Ks _1 m“ 1 . 
Taking a typical meridional temperature gradient of O.OIKkm -1 , we find that the zonal 
mean atmospheric cloud effect is equivalent to about 5% of the temperature gradient if 
clouds persist 4 to 5 days. It is, therefore, not negligible, although the meridional gradient 
of cloud effect is a small fraction of the rate of generation of mean zonal available potential 
energy. 

The zonal mean energy in the atmosphere changes with time because of the meridional 
energy flux gradient in the atmosphere, TOA zonally averaged net irradiance , and 
enthalpy flux from the surface, 


dE 

dt 


h dcosdFe tm 

' ^ rad 


R cos 0 dO 


( 8 ) 


where E is the zonal mean energy in a vertical atmospheric column, 0 is the latitude, Fq 
is the vertically averaged meridional component of the energy flux, R is the radius of the 
earth, and h is the height of the atmosphere. In a steady state condition, the radiative and 
enthalpy fluxes are balanced with meridional energy flux divergence in the atmosphere 

h dcosOFe = atm 
RcosO 80 rad H 
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When F + Fjj is integrated along latitude, we obtain the rate of energy crossing a 
vertical latitudinal plane in the atmosphere (e.g. Zhang and Rossow, 1997; Trenberth and 
Caron, 2001; Trenberth and Stepaniak 2003). As discussed in the previous section, we 
assume that the zonal mean cloud effect on + Fh can be separated from those under 
clear-sky conditions such that F^™' + Fh = where = (F^'j 1 + F H ) c id 

and = (F^™ + Fh) c W- We then estimate the zonal mean cloud effect on the vertical 

mean meridional energy flux F c by setting 

vatm h 9 cos 6F C 

Fm ~~r co S « ae ' (10) 

Then, (8) is 


dE 


h c) cos 0(F e + F c ) tm 

f-^clr > l 11 ) 


where 


dt R cos 9 d6 

h d cos 6 F c - g | j 10 meridional convergence of the energy flux that is equivalent 


R cos 6 80 

to the zonal mean cloud effect. Equation 11 can be considered as the energy equation 
for a clear-sky atmospheric column. Under clear-sky conditions, Ay j s no |, zero because 
the temperature profile needs to adjust from that under all-sky conditions. The global 
mean rate of the energy change in (11) is negative because the global mean atmospheric 
cloud radiative effect is slightly positive (~ 2Wm -2 ). To reach a steady state, therefore, 
the atmosphere needs to adjust the temperature profile corresponding to the meridional 
clear-sky net irradiance distribution. 

To obtain F c , we substitute (9) in (11), 

h d cos 9F C 


rp atm | ~ ^ 

— * eld, + "T7T- 


dE 

~dt 


( 12 ) 


R cos 9 d9 

Here, we are not interested in the meridional profile after the system reaches a steady 


state under clear-sky conditions. We are interested in estimating the energy transport 
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that is equivalent to the zonal mean cloud effect under all-sky conditions. Therefore, (12) 
is integrated over latitude and longitude by setting that — ^ equal to the seasonal global 
mean of F£{™ for all latitudes so that F c at poles is zero. The resulting 2ir RhF c cos 6 is the 
rate of the energy crossing a vertical latitudinal plane in the atmosphere that is equivalent 
to the zonal mean cloud effect. The underlying assumption is that the local rate of the 
energy change is uniform when clouds are removed. 

Figure 10 shows the rate of meridional energy transport in the atmosphere 
2ir RhF c cos 0 separated by seasons. Positive values indicate northward transport and 
negative values indicate southward transport. Two lines that provide the uncertainty 
envelope are shown in each plot. The solid line is the total cloud effect integrated along 
latitude and the dash-dot line considers only the radiative effect. As the effect of 
radiation on the surface enthalpy flux increases, 2 nRhF c cos 9 approaches to the line that 
includes the surface enthalpy change in the cloud effect. The energy transport has a 
maximum (minimum) in the northern hemisphere (southern hemisphere) midlatitude and 
it is mostly positive (negative) in the midlatitude to polar regions, especially in the 
fall and winter hemisphere. This estimate neglects dynamical feedbacks such as the 
enthalpy flux change due to wind speed change suggested by Randall et al. (1989). 
When the cloud indirect effect on the surface enthalpy flux is included, however, the 
study by Randall et al. (1989) suggests that the cloud effect on the meridional enthalpy 
transport might be larger in the tropics. These two estimates (dash-dot and solid lines in 
Figure 10), therefore, can be considered as an envelope of the rate of the meridional 
atmospheric energy transport equivalent to the zonal mean direct cloud effect. This 
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result indicates that the meridional gradient of cloud effect increases the rate of poleward 
energy transport from midlatitude to polar regions by dynamics in the atmosphere. 

As mentioned earlier, the atmospheric temperature under a clear-sky condition can be 
different from that under the all sky condition. A clear-sky atmosphere in midlatitude and 
polar regions can be warmer so that the net atmospheric irradiance under all-sky conditions 
may not be much different from the clear-sky value. This is unlikely for two reasons. First, 
a clear-sky atmosphere needs to be warmer by more than 10 K to compensate the cloud 
effect of about 30Wm -2 regionally. Second, the vertical distribution of cloud radiative 
effect works favorably for dynamics to transport energy poleward for the following reason. 
In a case where clouds move with winds, we apply the thermodynamic energy equation to 
a column moving poleward 


DE 

~Di 


v £p 

Dt 


T Fradi 


(13) 


where V is the volume of the column. When we neglect the pressure change with time and 
integrate this over time we obtain the energy transported by the column 


1 f yi 

[Eih) - E(t 0 )] = - F c i d (y ) + F clr (y)dy, (14) 

v Jyo 

where v is the meridional velocity of the column and the column is at y () at. t — to and 
at yi at t = t\. The column is heated or cooled by radiation F c id + F c i r at y. If clouds 
cool the column when they move poleward with the column, E(t i) is lower than that in 
a column without clouds so that the cloud reduces the energy transport. Therefore, the 
effect of clouds on the energy transport depends on the movement of clouds relative to 
wind. Figure 7 suggests that, however, the atmospheric cooling effect is caused by low-level 
clouds. Low-level clouds tend to form locally and remain relatively stationary. Because a 


21 



stronger cooling in polar regions is caused by low-level clouds, which subsequently causes 
the meridional temperature gradient, the rate of meridional energy transport to polar 
regions increases. High- and middle-level clouds, which tend to be advected by wind, 
have a warming effect to the atmosphere (Figure 7). When high- and middle-level clouds 
are advected poleward, the covariance of the column integrated temperature deviation 
due to clouds and the meridional component of velocity is positive. Because mean zonal 
available potential energy is converted to eddy available potential energy through poleward 
transport of warm air (Lorenz, 1955; Peixoto and Oort 1992, p377), this indicates that 
the rate of conversion of mean to eddy available potential energy is also increased when 
mid and high-level clouds are advected. Therefore, the meridional and vertical gradients 
of cloud effect increase the rate of generation of mean zonal available potential energy, rate 
of conversion of mean to eddy available potential energy, and rate of meridional energy 
transport from midlatitude to polar region through the atmosphere. 

5. Conclusions 

We estimated the zonal mean radiative cloud effect to the atmosphere using three years of 
CERES data. Although it is mostly positive for all four seasons, the zonal mean shortwave 
atmospheric cloud effect is small. A large effect occurs near midlatitudes in spring and 
summer presumably due to a large cloud fraction. Clouds, therefore, increase the shortwave 
absorption in the atmosphere. The net zonal mean atmospheric cloud radiative effect is, 
however, dominated by the longwave effect. The zonal mean longwave effect is a warming 
in the tropics and decreases approximately linearly with latitude to a cooling effect in 
polar regions. The meridional dependent atmospheric cloud radiative effect is a result 
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of decreasing the TOA cloud effect and increasing the surface cloud effect with latitude. 
This meridional dependence of the cloud effects is caused by the latitudinal difference of 
the mean cloud top height and water vapor effective emission height determined from the 
surface and the latitudinal variation in the cloud fraction. 

We consider the upper limit of the cloud effect of the surface enthalpy flux by the 
assumption that the cloud effect of the surface enthalpy flux is proportional to the net 
irradiance at the surface. The uncertainty envelope of the atmospheric cloud effect is 
between the upper limit and lower limit that is simply given by no cloud effect on the 
surface enthalpy flux. While the uncertainty in the cloud effect on the surface enthalpy 
flux is large in the tropics, the meridional gradient between midlatitude and polar regions 
exists even when uncertainties in the cloud effect on the surface enthalpy flux and in the 
modeled irradiances are taken into account. 

The atmospheric cooling effect of clouds in midlatitude and polar regions is caused by 
low-level clouds. Frequently occurring low-level clouds and their stronger cooling effect in 
polar regions lead to the meridional gradient of cloud effect between midlatitude and polar 
regions. The meridional gradient of cloud effect increases the rate of generation of mean 
zonal available potential energy. Because the cooling effect over polar regions is caused 
by relatively stationary low-level clouds, the meridional gradient of cloud effect increases 
the rate of meridional energy transport. Middle- and high-level clouds, which tend to be 
advected by winds have warming effects to the atmosphere. This leads to increase the rate 
of conversion of mean to eddy available potential energy. Clouds in polar regions then 
warm the surface except in the Arctic in summer. 
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Because clouds form as a consequence of dynamics, these results suggest a possible 
feedback process (e.g. Weaver 2003) among cloud type, meridional energy transport, and 
surface temperature in polar regions. If we consider that about 50% of the energy emitted 
to space by polar regions is provided from midlatitudes (Kato et al. 2006), understanding 
the link between energy transport and the zonal mean cloud effect is as equally important 
as understanding local cloud feedback processes in polar regions. 
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Appendix A: Daily Mean Irradiance Estimate 

To convert instantaneous irradiances to a daily value, we need to obtain the TOA albedo, 
transmission, and surface albedo as a function of solar zenith angle. We sort and average 
instantaneous TOA albedos estimated from CERES radiance as a function of solar zenith 


angle and scene type using data taken from the TRMM satellite from January 1998 through 
August 1998 and March 2000. Note that the effect of a spherical earth that leads to a lion- 
negligible TOA shortwave irradiance over regions where the solar zenith angle is greater 
than 90° is included based on Kato and Loeb (2003). Similarly, we sort and average 
transmissions and surface albedos as a function of solar zenith angle and scene type using 
one- year of Terra CRS from March 2000 through February 2001. The increment of the 
solar zenith angle bin is 10° for the TOA albedo and 5° for the surface transmission and 
albedo. The number of scene types that contain no snow and sea ice is 590 (Loeb et al. 
2003). The number of snow and sea ice scenes is 60 (Kato and Loeb 2005). 

The instantaneous irradiance F^^ A is scaled by the ratio of the daily mean TOA 
albedo A TOA (x ) to the instantaneous TOA albedo Af OA (9o, x), 

Fj° A (x) = Fl° w A (0 o ,x) [ > (al) 

Af UA (9 0 ,x) 

where Fj^ A is the daily mean TOA irradiances, 9q is the solar zenith angle, and x is the 
scene type over a CERES footprint that is estimated from MODIS radiances. The diurnal 
change of solar zenith angle over the course of a day for the location of the radiance 
measurement is used to compute A , OA (x). Similarly, the daily mean downward and 
upward shortwave surface irradiances are computed by 


F sw^dn ( X ) = F f,dn( 0 O,x) 


Ti(9 0 , x y 
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and 




A(x) 

Ai(0 0 ,x)’ 


(a3) 


where T{x) and A(x) are transmittance and surface albedo averaged over the course of a 
day for the scene type x, and Ti(9o,x) and Ai(9o,x) are the instantaneous transmittance 
and surface albedo at the CERES overpass time, respectively. 


Appendix B: Zonal Mean Atmospheric Cloud Radiative Effect 

Figure 6 shows that the zonal mean atmospheric cloud radiative effect, which is dominated 
by the longwave effect, changes almost linearly with latitude. In addition, the meridional 
gradient is larger than the seasonal variations. In this appendix, we explain the reason for 
this simple latitudinal dependence of the zonal mean atmospheric cloud radiative effect. 
Because the atmospheric cloud radiative effect is dominated by longwave, we only consider 
the longwave effect here. We utilize a plane parallel atmosphere with the temperature T a , 
and emissivity e. The temperature of the underlying surface is T s and the emissivity is 
assumed to be unity. The top-of-atmosphere longwave net irradiance F^ OA and surface 
longwave net irradiance F'lJf is 

F i ™- 4 = -[(l-£) C rT 5 4 + Et rr o 4 ], (04) 


and 


Ft = ^r 4 


aT;. 


(a5) 


If we let T a — T s + AT, then 


7 toa 

Iw 


-aT s 4 (l + 


4eAT, 
T, ' 


(a6) 
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and 


F sfc « 

r lw 


-vTsi 1 ~ e - 


4eAT 
T J 

-*■ s 


(a7) 


The longwave cloud radiative effect at the top of the atmosphere F^ OA — Fu,?cir an d the 
surface F^ c - F^ clr is 


rpTOA _ rpTOA 
Iw ^ lw,clr 


= -4crT s 3 (eAT - e cir AT cir )/ c 


(a8) 


and 


ps/c _ ps/c 
± Iw ± lw,clr 


<rr 5 4 [€(l+4 



— £clr{ 1 + 4 


AT e | r 

T s 


)]fc, 


(o9) 


where the subscript c/r indicates the clear condition and f c indicates cloud fraction. Be- 
cause clouds reduce effective temperature when viewed from TOA, eA T — e c i r AT c i r is 
negative. When viewed from the surface, clouds increase the emissivity of the atmosphere 
and effective temperature, the right side of (a9) is positive. 

If we assume that the difference between the effective temperature and surface tem- 
perature AT is 


(IT 

AT = — z top , 

_ dT " 

AA-L clr Ztop,WV 5 


and the emissivity of clear-sky atmosphere is 


(alO) 

(all) 


e c ir = 1 - e r ’ 


(al2) 


where z is the altitude, zt op is the cloud top height, zt op ,wv is the effective height of water 
vapor that gives effective temperature of clear-sky when viewed from TOA, and r wv is the 
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optical thickness of water vapor. We also assume that the emissivity of clouds is unity, 
e = 1. In addition, we assume 


A T _dT 

" base: 

clz 


AT - <U L 

clr 7 ^base.wv 5 

az 


(ol3) 

(o!4) 


and 


^ clr 


1 — e 


7~WV 

5 


(o!5) 


where Zbase is the cloud base height and Zbase, wv is the effective height of water vapor that 
gives effective temperature of clear-sky when viewed from the surface. 

Substituting (alO) through (al5) into (a8) and (a9) and assuming that the column 
water vapor amount is sufficient to give r wv 3> 1, we then obtain 

fL oa - FL°cl = - z top , wv )f c , (al6) 

and 

dT 

F lw° - F hl%r = ° T S -fa ( Zbase ~ Zbase, wv)fc («17) 

Figure 5 shows that the meridional gradient of atmospheric cloud radiative effect in 
tropics and midlatitudes is caused by both TOA and surface effects and that over polar 
regions is predominately caused by the TOA effect. Figures al and 3 indicate that zt op 
decreases with latitude between 20N and 20S probably more than z top ,wv does. The cloud 
fraction f c also decreases with latitude between 20N and 20S. Therefore, the zonal mean 
TOA cloud radiative effect decreases with latitude in this region. Because the column water 
vapor amount decreases with latitude, Zbase, wv increases while Zbase is relatively constant 
with latitude. As a consequence, the zonal mean cloud radiative effect at the surface 
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increases with latitude. Figure al also indicates that the cloud height decreases with 
latitude in polar regions. Therefore the zonal mean TOA cloud radiative effect decreases 
with latitude in polar regions. Because the zonal mean atmospheric cloud radiative effect 
is (al7) subtracted from (al6), 

( F L OA - F hu°cir )-( F lw C ~ F lw%r) = JZ [i z top ~Ztop,wv) + {z base ~ Zbase, wv)] fc- (0I8) 


The zonal mean atmospheric cloud radiative effect is, therefore, a function of the difference 
between the mean cloud height and water vapor effective emission height multiplied by the 
cloud fraction. In the tropics, Zt op is higher than Zt op ,wv and Zbase is nearly equal to 
Zbase, wv Therefore, the zonal mean atmospheric cloud radiative effect is positive in the 
tropics. In polar regions, the difference between Zt op and Zt op ,wv is small while Zbase is 
less than Zbase, wv Therefore, the zonal mean atmospheric cloud radiative effect is negative 
in polar regions. If we assume the difference between Zbase and Zt op ,wv is small compared 
with the difference between z top and Zb as e,wv> then 


(Fiw 


TOA 


7 TOA 

lw,clr 


) - ( F i 


sfc 


7 <sfc 
Iw.clr 


) = 


, rT ~. q dT 
- iaT ° Tz {z ‘ 


op 


Zbase, wv )fc 


(o!9) 


which indicates that the zonal mean atmospheric cloud radiative effect is, to a first order 
approximation, a function of difference between the mean cloud top height and water vapor 
effective emission height determined from the surface longwave irradiance multiplied by 
the cloud fraction. While this does not explain the reason for a simple linear relation with 
latitude, it tells that the sign of the zonal mean atmospheric cloud radiative effect changes 
going from the tropics to polar regions. 
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Figure 1 Seasonal and zonal mean net shortwave irradiance a), net longwave irradiance b), 
and net shortwave and longwave irradiance c) of the atmosphere under all-sky conditions. 
The atmospheric shortwave absorptance d) is computed from a) divided by zonal mean 
solar constant. The net longwave irradiance and net shortwave plus longwave irradiance 
of the atmosphere under clear-sky conditions are shown by e) and f), respectively. 

Figure 2 Comparison of monthly mean downward shortwave (top), longwave (middle) and 
net atmospheric (bottom) irradiances for Manus (TWP), Southern Great Plains (SGP), 
and Barrow, AK (NSA) sites. Open circles and closed squares indicate modeled irradiances 
and observations derived from March 2000 through Feb. 2003, respectively. The error bars 
indicate maximum and minimum observed values during the 5-year period (March 2000 
through Feb. 2005) for the surface down shortwave and longwave and 4-year period (March 
2000 through Feb. 2004) for the atmospheric net. 

Figure 3 a) Seasonal zonal mean cloud fraction derived from MODIS measurements by the 
CERES cloud algorithm (Minnis et al. 2003). Data from March 2000 through Feb. 2004 
are averaged, b) Monthly zonal mean cloud fraction derived from GLAS 1064nm laser 
October 2003 data for all clouds (dotted line), those of which optical thickness greater 
than 0.3 (solid line) and optical thickness greater than 1 (dash-dot line). The cloud optical 
thickness is derived from GLAS 532 nm laser data. The cloud fraction derived from 
MODIS radiances for the same month is shown by the dashed line, c) Difference between 
the MODIS derived cloud fraction and the GLAS derived cloud fraction. The fraction of 
clouds of which optical thickness is less than 0.3 is excluded from the GLAS derived cloud 
fraction. 
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Figure 4 Seasonal and zonal mean cloud shortwave radiative effect at TOA a), at the 
surface b), and to the atmosphere c). 

Figure 5 Seasonal and zonal mean cloud longwave radiative effect at TOA a), at the 
surface b), and to the atmosphere c). 

Figure 6 Seasonal and zonal mean cloud radiative effect (shortwave plus longwave) at 
TOA a), at the surface b), and to the atmosphere c). 

Figure 7 Contour of the daily mean cloud shortwave plus longwave effect at TOA (top 
row) to the atmosphere (second row) and to the surface (third row) for the tropics (30N to 
30S), northern hemisphere mid-latitude (30N to 60N), southern hemisphere mid latitude 
(30S to 60N), the Arctic (60N to 90N) and the Antarctic (60S to 90S) as a function of 
the cloud optical thickness (r) and cloud top height in pressure coordinate estimated from 
July 2002 data. Only single-layer clouds are used. Contours in the forth row indicate 
the logarithm (base 10) of the 2D normalized histogram of cloud occurrence. Daily mean 
irradiances are computed by the method discussed in Appendix A except that daytime 
and nighttime longwave irradiances are weighted by number of samples. 

Figure 8 Normalized 2D histogram sorted by net surface irradiance (shortwave+longwave) 
anomaly divided by monthly mean value and surface enthalpy (latent heat and sensible 
heat) flux anomaly divided by monthly mean value. The latent heat and sensible heat 
fluxes are from Hamburg Ocean Atmosphere Parameters and fluxes from Satellite data 
(HOAPS, Grassl et al. 2000). The anomalies are defined as the deviation of the zonal 
mean values from the averaged value over the four-year period from March 2000 through 
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Dec. 2004. One month of data (March 2000) between 20°N and 20° S were used for the 
plot. The solid line indicates linear regression fit. The the slope of the regression line is 
correlation coefficient is 0.16. 

Figure 9 Seasonal and zonal mean cloud effect to the atmosphere. The effect includes that 
on the shortwave and longwave irradiances, as well as on the surface latent and sensible 
heat fluxes (enthalpy flux). 

Figure 10 Seasonal and zonal mean atmospheric cloud effect converted to the rate of 
meridional energy transport by the atmosphere. The solid line indicate the equivalent 
meridional energy transport including the cloud effect on the net irradiance, surface latent 
heat, and sensible heat fluxes. The dash-dot line indicates the equivalent meridional en- 
thalpy transport including only the cloud effect on the net irradiance. The positive and 
negative value indicates, respectively, northward and southward transport. 

Figure al Monthly and zonal mean cloud top pressure derived from daytime July 2002 
MODIS data by the CERES cloud algorithm (Minnis et al. 2003). Only single layer clouds 
are considered. 
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Table 1: Observed and Modeled Annual Mean Irradiance 


Observation 

Wm“ 2 

Model 

Wm“ 2 

Difference 

Wm“ 2 

Relative 

Difference 

Manus (Tropics) 

SW TOA Down 

417 

417 

- 

- 

SW TOA UP 

110 

111 

1.1 

0.010 

SW Surface Down 

208 

228 

20.5 

0.099 

SW Surface Up 

69 

17 

- 

- 

SW Atm. Net 

117 1 

95 

-20.7 

-0.177 

LW TOA UP 

219 

220 

1.3 

0.006 

LW Surface Down 

420 

419 

-0.6 

-0.001 

LW Surface Up 

469 

476 

7.8 

0.017 

LW Atm. Net 

-170 

-163 

7.1 

0.041 

SW+LW Atm. Net 

-53 1 

-68 

-14.6 

-0.274 

SGP (Midlatitude) 

SW TOA Down 

343 

343 

- 

- 

SW TOA UP 

103 

104 

1.0 

0.010 

SW Surface Down 

187 

199 

12.6 

0.068 

SW Surface Up 

39 

39 

-0.5 

0.008 

SW Atm. Net 

92 

78 

-14.1 

-0.154 

LW TOA UP 

243 

243 

0.2 

0.008 

LW Surface Down 

337 

326 

-10.8 

-0.032 

LW Surface Up 

395 

402 

6.6 

0.017 

LW Atm. Net 

-185 

-168 

17.2 

0.093 

SW+LW Atm. Net 

-93 

-90 

3.1 

0.033 

Barrow (Arctic) 
SW TOA Down 

196 

196 



SW TOA UP 

95 

98 

3.0 

0.032 

SW Surface Down 

96 

102 

6.4 

0.067 

SW Surface Up 

48 

49 

1.3 

0.027 

SW Atm. Net 

53 

45 

-8.1 

-0.153 

LW TOA UP 

203 

205 

1.8 

0.009 

LW Surface Down 

238 

238 

-0.3 

0.001 

LW Surface Up 

271 

271 

0.0 

0.000 

LW Atm. Net 

-171 

-172 

-1.6 

-0.009 

SW+LW Atm. Net 

-118 

-127 

-9.7 

-0.082 


1 


Modeled shortwave upward surface irradiance is used for the computations. 
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Figure 1 Seasonal and zonal mean net shortwave irradiance a), net longwave irra- 
diance b), and net shortwave and longwave irradiance c) of the atmosphere under all-sky 
conditions. The atmospheric shortwave absorptance d) is computed from a) divided by 
zonal mean solar constant. The net longwave irradiance and net shortwave plus longwave 
irradiance of the atmosphere under clear-sky conditions are shown by e) and f), respec- 
tively. 
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Figure 2 Comparison of monthly mean downward shortwave (top), longwave (mid- 
dle) and net atmospheric (bottom) irradiances for Manus (TWP), Southern Great Plains 
(SGP), and Barrow, AK (NSA) sites. Open circles and closed squares indicate modeled 
irradiances and observations derived from March 2000 through Feb. 2003, respectively. 
The error bars indicate maximum and minimum observed values during the 5-year period 
(March 2000 through Feb. 2005) for the surface down shortwave and longwave and 4-year 
period (March 2000 through Feb. 2004) for the atmospheric net. 
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Figure 3 a) Seasonal zonal mean cloud fraction derived from MODIS measurements 

by the CERES cloud algorithm (Minnis et al. 2003). Data from March 2000 through Feb. 

2004 are averaged, b) Monthly zonal mean cloud fraction derived from GLAS 1064nm 

laser October 2003 data for all clouds (dotted line), those of which optical thickness greater 

than 0.3 (solid line) and optical thickness greater than 1 (dash-dot line). The cloud optical 

thickness is derived from GLAS 532 nm laser data. The cloud fraction derived from 

MODIS radiances for the same month is shown by the dashed line, c) Difference between 

the MODIS derived cloud fraction and the GLAS derived cloud fraction. The fraction of 

clouds of which optical thickness is less than 0.3 is excluded from the GLAS derived cloud 
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Figure 4 Seasonal and zonal mean cloud shortwave radiative effect at TOA a), at the 


surface b), and to the atmosphere c). 
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Figure 5 Seasonal and zonal mean cloud longwave radiative effect at TOA a), at the 


surface b), and to the atmosphere c). 
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Figure 6 Seasonal and zonal mean cloud radiative effect (shortwave plus longwave) 
at TOA a), at the surface b), and to the atmosphere c). 
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Figure 7 Contour of the daily mean cloud shortwave plus longwave effect at TOA 
(top row) to the atmosphere (second row) and to the surface (third row) for the tropics 
(30N to 30S), northern hemisphere mid-latitude (30N to 60N), southern hemisphere mid 
latitude (30S to 60N), the Arctic (60N to 90N) and the Antarctic (60S to 90S) as a function 
of the cloud optical thickness (r) and cloud top height in pressure coordinate estimated 
from July 2002 data. Only single-layer clouds are used. Contours in the forth row indicate 
the logarithm (base 10) of the 2D normalized histogram of cloud occurrence. Daily mean 
irradiances are computed by the method discussed in Appendix A except that daytime 
and nighttime longwave irradiances are weighted by number of samples. 
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Surface Net Irradiance Relative Anomaly 

Figure 8 Normalized 2D histogram sorted by net surface irradiance 
(shortwave+longwave) anomaly divided by monthly mean value and surface enthalpy 
(latent heat and sensible heat) flux anomaly divided by monthly mean value. The latent 
heat and sensible heat fluxes are from Hamburg Ocean Atmosphere Parameters and 
fluxes from Satellite data (HOAPS, Grassl et al. 2000). The anomalies are defined as the 
deviation of the zonal mean values from the averaged value over the four-year period 
from March 2000 through Dec. 2004. One month of data (March 2000) between 20°N 
and 20°S were used for the plot. The solid line indicates linear regression fit. The the 
slope of the regression line is correlation coefficient is 0.16. 
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Figure 9 Seasonal and zonal mean cloud effect to the atmosphere. The effect includes 


that on the shortwave and longwave irradiances, as well as on the surface latent and sensible 


heat fluxes (enthalpy flux). 
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Figure 10 Seasonal and zonal mean atmospheric cloud effect converted to the rate 
of meridional energy transport by the atmosphere. The solid line indicate the equivalent 
meridional energy transport including the cloud effect on the net irradiance, surface latent 
heat, and sensible heat fluxes. The dash-dot line indicates the equivalent meridional en- 
thalpy transport including only the cloud effect on the net irradiance. The positive and 
negative value indicates, respectively, northward and southward transport. 
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Figure al Monthly and zonal mean cloud top pressure derived from daytime July 
2002 MODIS data by the CERES cloud algorithm (Minnis et al. 2003). Only single layer 
clouds are considered. 
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